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Abstract. We study the single-species diffusion-annihilation process with a 
time-dependent reaction rate, A(t) = Aq Scaling arguments show that there is 
a critical value of the decay exponent LJc{<i) separating a reaction-limited regime 
for w > uJc from a diffusion-limited regime for uj < uJc- The particle density 
displays a mean-field, oj-dependent, decay when the process is reaction limited 
whereas it behaves as for a constant reaction rate when the process is diffusion 
limited. These results are confirmed by Monte Carlo simulations. They allow 
us to discuss the scaling behaviour of coupled diffusion-annihilation processes in 
terms of effective time-dependent reaction rates. 



PACS numbers: 82.20.-w, 05.40.-a, 05.70.Ln 

1. Introduction 

Time-dependent effective reaction rates can be introduced into mean-field kinetic 
equations in order to simulate the effect of concentration fluctuations below the upper 
critical dimension where deviations from the standard mean-field behaviour are 
observed f . In this work we take the opposite point of view and look how far 
the fluctuations in a reaction-diffusion process can be affected by a time-dependent 
reaction rate. 

We study the single-species diffusion-annihilation process 
A0 0A (diffusion) 

A A ^ (annihilation), (I) 
where nearest-neighbour particles annihilate with a reaction rate 

m = (2) 

decaying as a power of the time. The diffusion rate associated with the exchange of a 
particle (A) with a vacancy (0) is equal to I, which fixes the time scale. 

When the reaction rate is constant {u; — 0), the kinetic equation p = — Aop^ leads 
to the mean-field asymptotic behaviour p{t) ~ (Aoi)~^. Simple scaling arguments [2E1 
indicate that, due to concentration fluctuations, the decay is slower in low dimensions. 
At time t, for d < 2 where the exploration is dense, a surviving particle has swept out 



Reaction- diffusion with a time- dependent reaction rate 



2 



a region with a linear size given by the diffusion length V Dt. Therefore, the volume 
per particle grows as (Dt)'^^^ and the particle density decays algebraically as 

p(t)=C7t-", d<dc-2. (3) 

Thus mean-field theory gives the correct scaling behaviour above the upper critical 
dimension d^ — 2. 

This picture is confirmed by exact results in one dimension (fd) [4-6], rigourous 
bounds [7| and renormalization group results |S1 IHI showing that the amplitude C is 
universal. As usual, logarithmic corrections occur at dc where the particle density 
decays as 

, , I hit , , „ , ,^ 

Pit) = ^ , d = d, = 2, 4 
Stt Dt 

where D is the diffusion constant. 

We consider the influence of the decay exponent uj, the amplitude Aq and the 
dimension d of the system on the scaling behaviour of the particle density p{t). In 
section 2, we give scaling arguments showing that the mean-field behaviour should 
be recovered when the decay exponent is greater than a critical value ujc{d) when 
d < dc and calculate the time evolution of the particle density in mean-field theory. In 
section 3 these results are confronted with Monte Carlo simulations data in dimension 
d = 1 to 3. In section 4 we show how time-dependent reaction rates are effectively 
realized in the case of coupled reactions. Our results are summarized in section 5. 



2. Scaling considerations and mean-field theory 

In a continuum description in d dimensions, the spacetime evolution of the particle 
density field, p(r, t), is governed by a Langevin equation containing diffusion, reaction 
and noise terms [10-12] 

dtpir, t) = p^r, t) +DV^p{r, t) + C(r, t) . (5) 

The noise term C(r, t) accounts for the fluctuations of the particle density at position 
r at time t. 

Let us consider the behaviour under rescaling of the reaction term 

9tp(r,i) [reaction = -^/0^(r, i) (6) 

at the stable fixed point of the system with a constant reaction rate, when the 
fluctuations are relevant, i.e., below the upper critical dimension dc- 

Under a change of the length scale L' ~ L/b, the particle density fleld and its 
space average p(t) scale with the same dimension Xp and the scaling dimension of t is 
the dynamical exponent z = 2, so that 

[pir,t)r^b^^pir,t), t'^-L. (7) 

For the average particle density, one obtains 

Taking 6 = t^/"^ leads to the power-law decay 

p{t)^p{l)t-\ a = ^. (9) 
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Using the transformations Q in ||HJl, we obtain 

9t'[/9(r,i)]' Inaction = &^''+''9t/9(r, i) Inaction 

= -¥-'^+^^'-^ p^ir,t) (10) 
so that the reaction-rate amplitude transforms as 

where the last expression follows from 0. 

When (jj is smaller than the critical value Uc, given by 

Wc = l- a = l- ^, d<dc, (12) 

according to (PJ, Xq increases under rescaling, i.e., the process is diffusion limited. 
The concentration fluctuations are relevant and the critical behaviour is governed by 
the same fixed point as for a constant reaction rate, for which a = d/2. 

When Lo > LOc, the reaction-rate amplitude decreases and the process is reaction 
limited. The concentration fluctuations are suppressed by diffusion and the system 
should display a mean- field behaviour governed by a fixed line, parametrized by to. 

When w = Wc, Ao is a marginal variable and logarithmic corrections to the mean- 
field behaviour are expected. 

The behaviour of p{t) follows from the mean-field rate equation 

p{t) = -\{t)p'it) (13) 
With the initial condition p{to) = pQ, the average density reads 

p{t) po 1-uj 
^ =- + Xoln(f) icj = l) (14) 



p{t) Po \to 
In the asymptotic regime, t S> ioi one obtains 

p{t) ~ (^<i) (15a) 

Aq 

Pit) ^ 3^ = (15b) 

^-(t^-i) 

p{t) ~ Poo+pIo Ao > 1) (15c) 

U! — 1 

When Lo > 1, the particle density decays algebraically to a non- vanishing asymptotic 
value. Poo, behaving as 



1 1 . t 



— = - + Xo^ (cc>l) 

Poo Po LU - I 

Poo =i (t^ ^ 1 + ) (16) 

When d > dc, ujc = according to H12|) . For positive values of w, the annihilation 
process is reaction limited and the mean-field behaviour in (15a)-(15c) and (fTB|l 
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applies. For negative values of w the reaction rate, which increases under rescaling, 
no longer controls the annihilation process. Thus one expect a scaling behaviour 
independent of w, the same as for a; = with a ~ I and logarithmic corrections given 
by ^ at dc when w < 0. 

3. Monte Carlo simulations 

Numerical simulations of the diffusion-annihilation process have been performed in 
order to check the results of the last section for the asymptotic behaviour of the 
particle density p(t) in dimensions d = 1 to 3. 

3.1. Algorithm 

We work on hypercubic lattices with L"^ = 10^ sites and periodic boundary conditions 
in all directions. In order to avoid finite-size effects the simulations are stopped when 
the diffusion length reaches some fraction (1/4 to 1/10) of the size L of the system. 
The particle density p{t) is averaged over 5 to 20 samples. 

At io = 1 the L'^ sites are independently occupied by a particle with probability 
Pq. At time t, one of the N{t) surviving particles is randomly selected and a jump 
towards one of the 2d neighbouring sites is attempted, with the same probability for all 
the sites. When the target site is empty, the jump is accepted. When it is occupied, 
either the two particles annihilate with probability X{t) or they keep their original 
location with probability 1 — A(t), i.e., multiple occupancy of a site is not allowed. 
Finally, the time is incremented by 1/N{t) in all the cases and the process is repeated. 

When cj < 0, the algorithm has to be modified since the reaction rate may increase 
beyond 1. In this case, we let two particles annihilate with probability \{t)T{t) = 1. 
We take Aq = 1 such that \{t) > 1 and T{t) < 1. The diffusion jumps are attempted 
with probability r(t) < 1 and the time is incremented by T{t)/N{t) at each Monte 
carlo step. 

3.2. Numerical results in Id 

The influence of cu on the scaling behaviour of p{t) is shown in figure^ The asymtotic 
slope, —a, is the same as for a constant reaction rate (—1/2, indicated by a dashed 
line) as long as a; < Wc = 1/2. The amplitude of p(t) is independent of w below ujc 
and it is universal as for lu — |^ : it depends neither on the initial density pq nor on 
the reaction-rate amplitude Aq, the process being diffusion limited. 

The form of the logarithmic correction to the algebraic decay, expected at Uc in 
relation with the marginality of Aq , could not be extracted from our Monte Carlo data. 

Above ujc, a decreases continuously to when uj increases to 1. As shown in 
figure[21the amplitude of p{t) is no longer universal: it remains independent of po but, 
the process being now reaction limited, it depends on Aq. The inset shows that the 
amplitude actually varies as 1/Ao, in agreement with the mean-field result in (|15b .). 

The scaling behaviour of p{t) at w = 1 is illustrated in figure |3| The product 
plii{t) tends to a constant value at long times, in agreement with mean-field theory 
in equation IjlSb ). 

The variation of a with uj is illustrated in figure 01 The exponents were deduced 
from the extrapolation of two-point approximants for the slope of ln(/5) versus hi{t) 
using the Burlisch-Stoer (BS) algorithm \l'6l I14| . These data were obtained with 
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ln(t) 

Figure 1. Time dependence of tlie particle density in Id for different values of tlie 
decay exponent lo. The initial density is p = 0.3 and the reaction-rate amplitude 
Ao = 1. The asymptotic slope, —a, varies with lo in the reaction-limited mean- 
field regime, above the critical value uic = 1/2. Below lUc, in the diffusion-limited 
regime, it remains constant and equal to —1/2 (dashed lines). The amplitude is 
universal below uic, the data for lo = 0.25, and —0.25 collapsing asymptotically 
on a single line. 




-6 I ■ • ■ • ' • ' • ' 1 
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Figure 2. Influence of po and Aq on the amplitude of the particle density p{t) 
in Id. At u> = 0.7, i.e., in the reaction-limited regime, the amplitude does not 
depend on the initial density po and it varies as 1/ Aq as shown in the inset where 
an asymptotical collapse is obtained for pAo. 



Po = 0.3 in order to accelerate the approach of the asymptotic regime (see figure |21 
for a comparison with po — I when Aq — 1). The value Aq = 1 was selected to spare 
computer time although smaller values can lead to better estimates of the exponent as 
illustrated for oj = 0.7 where the circles correspond to values of Aq — 1, 0.4, 0.2 and 0.1 
from top to bottom. The results are in good agreement with the expected behaviour. 
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Figure 3. Logarithmic scaling behaviour at oj = 1 for ci = 1 to 3. In the mean- 
field reaction-limited regime the particle density decays asymtotically as l/ln(p) 
when LU = I. The data were obtained with po = 1 and Aq = 1 
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Figure 4. Variation of the particle density exponent a with the reaction-rate 
exponent u! for ci = 1,2 and 3. The data, obtained with p = 0.3 and Aq = 1, 
are in overall agreement with the expected behaviour (solid lines). Below uJc (1/2 
when d = 1, when d = 2,3) the process is diffusion limited and q is a constant. 
Above oJc the process is reaction limited and a decays as 1 — a; in agreement with 
mean-field theory. 

3.3. Numerical results in 2d and 3d 

The log-log plots for the time evolution of the particle density for different values of 
Lv are shown in figure [S] for d — 2 and figure for d ^ 3. As above the data were 
obtained with an initial density po — 0.3 and a rection-rate amplitude Aq = 1. In both 
cases the critical value of the reaction-rate exponent is Uc = 0. a varies with lu in the 
reaction- limited regime, lu > lOc, and remains constant in the diffusion- limited regime, 

to < UJc. 
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Figure 5. Time dependence of the particle density in 2d with po = 0.3 and Aq = 
1. The asymptotic slope depends on lj in the reaction-limited regime ut > utc = 
and remains constant in the diffusion-limited regime (oj = 0, —0.25, —0.5). There 
is a deviation from the expected slope, —1 (dashed lines), due to the logarithmic 
correction occuring in the diffusion-limited regime at dc = 2. When p is divided 
by ln{t) (dotted lines) the correct slope is recovered. 




Figure 6. Time dependence of the particle density in 3d with po = 0.3 
and Ao = 1. The critical value of the reaction-rate exponent is lJc = 0, the 
same as in 2d. Since d > dc, there are no logarithmic corrections and the 
asymptotic slopes are close to —1 (dashed lines) in the diffusion-limited regime 
{uj = 0,-0.25,-0.5,-0.75). 



In 2d there is a systematic deviation from the slope —1 (indicated by a dashed 
line in figure |3J| when a; < 0. This may be traced to the logarithmic correction in 
equation occuring at the upper critical dimension dc = 2 in the diffusion-limited 
regime. When p is divided by \n{t) (dotted lines) the slopes are asymptotically close to 
the expected value —1 corresponding to the mean- field exponent a = I. The amplitude 
is probably universal, the same for all values of w < 0, although the evolution to the 
asymptotic behaviour becomes quite slow when approaching uJc — 0. 
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In d = 3 > dc the mean-field result of equation l^lbh ) with uj — 0, p{t) ^ t^^, 
applies in the whole diffusion-limited regime w < as shown in figure IHl where the 
slopes are close to the expected value, indicated by the dashed lines, without any 
correction. Here too we expect the amplitude to be universal when lu < but the 
evolution to the true asymptotic behaviour is even slower than in 2d when ujc is 
approached. 

The asymptotic logarithmic decay of the particle density at a; = 1 (see 
equation p5b ll is illustrated in figure|31 The variation of a with w, shown in figure^ is 
still in agreement with the expected behaviour. The exponents were obtained through 
a non-linear least-square fit, taking into account an effective correction-to-scaling 
exponent, since the statistical fluctuations were too high to use the BS algorithm. 

4. Time-dependent effective reaction rates 




Figure 7. Scaling behaviour of an asymmetric diffusion-annihilation process in 
2d where the annihilation of A particles is catalysed by B particles whereas the 
annihilation of B particles is free. The product pA lii^(t) tends to a constant value 
at long time. 



Time-dependent effective reaction rates with a power-law decay can be realized in 
the case of coupled reactions. Let us first consider a diffusion-annihilation process in 
2d where the annihilation of a pair of A particles is catalysed by B particles whereas 
B particles annihilate without any further condition: 

ABA — ^ B (catalysed annihilation) 

BB — ^ (simple annihilation) . (17) 

The effective reaction rate for AA annihilation is governed by the mean density 
of B particles, psit), if the concentration fluctuations are negligible, which is true in 
2d. Thus UJA — Ob — I for the single-species annihilation in 2d. The AA process 
being reaction limited with to a — 1, mean-field theory applies and a logarithmic decay 
is expected. Actually, one has to take into account the logarithmic correction in 
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Figure 8. Scaling beliaviour of tlie particle density for a 2d symmetric diffusion- 
annihilation process where the annihilation of two particles of one species is 
catalysed by a particle of the other species. The reaction rates are equal, 
r = Xa/^b = li and both densities decay as t~^/^. The density of B particles 
with initial density pg (0) is equal to the density of A particles with initial density 
Pa(0) = 1 — pb{0)- The asymptotic slope, —1/2, is indicated by a dashed line. 
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Figure 9. Scaling behaviour of the particle density for a 2d symmetric diffusion- 
annihilation process where the annihilation of two particles of one species is 
catalysed by a particle of the other species. Here the reaction rates are different 
with r = \a/^b = 1/2. The initial densities are such that Pa(0) + Pb(0) = 1 and 
equal to 0.7, 0.5, 0.3, from top to bottom for A and B particles. The expected 
asymptotic slopes, —1/3 for A particles and —2/3 for B particles, are indicated 
by dashed lines. The convergence is slow for the majority species. 



equation Q so that the effective reaction rate is XA{t) 
rate equation H13|) leads to the asymptotic behaviour 

2 1 

PA{t) 



AAoln(i)/*- The mean- field 



(18) 



Xao \n{t) ■ 

This process has been simulated on the square lattice with Xa — Xb — ^ and 
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Pyi(O) + pb(0) 1. The size of the system, the boundary conditions and the number 
of samples are the same as before. When a particle jump is attempted towards an 
occupied site, the two particles annihilate in the case of a i? pair. For a pair of A 
particles, the annihilation occurs only when a B particle is first neighbour of one of the 
A particles and second neighbour of the second. When the two particles are different 
they just keep their positions. 

The Monte Carlo results, shown in figure 7, confirm this scaling behaviour 
although the convergence is quite slow for small or intermediate initial densities of 
A particles. The second example concerns a system of two symmetrically coupled 
single-species diffusion-annihilation processes in 2d with 

ABA — ^ B (catalysed annihilation) 

BAB ^ 0A0 (catalysed annihilation) . (19) 

The density of one species controls the reaction rate of the other so that oja = ctB- 
Now, assuming that the system is in the reaction-limited regime, equation IjlSb l leads 
to a A = 1 — = 1 — Q!_B and 

aA + aB = l. (20) 

When the process is fully symmetric, i.e., when A^i — Xb, both species decay with 
the same exponent a a — cub — 1/2, a value consistent with the assumed reaction- 
limited mean-field behaviour since = in 2d. Actually, the fully symmetric process 
reduces to the 2>A — > A process which is known to decay as above its upper 

critical dimension dc = 1 |15| . 

The Monte Carlo simulations on the square lattice with r = A^/Ab = 1, 
/9^(0) -I- p_b(0) = 1 and the same rules as above for the catalysed annihilation of 
the two species are shown in figure 8. The scaling behaviour is in complete agreement 
with a decay exponent equal to 1/2. 

A deviation of the reaction-rate ratio r from 1 is sufficient to break the symmetry 
between the two species, leading to different decay exponents. This follows from the 
solution of the system of coupled mean-field rate equations 

f>A{t)^~XAPA{t)PB{t), pB{t) = -\Bpl{t)pA{t). (21) 

Multiplying the first equation by ps, the second by pA and adding, one obtains a 
differential equation for paPb leading to the asymptotic behaviour 

PA{t)pB{t) - t-^ (22) 

in agreement with H20|) . Dividing the first equation by pA and using (|22|l . one finally 
obtains 

PA{t)^t-'l^^+^\ pBit)^t-'/^'+-\ r = \A/\B. (23) 

The decay exponents continuously vary with the ratio r of the reaction rates but their 
sum remains equal to 1 as expected from the effective reaction-rate argument leading 
to (|20|l . Actually, it is easy to verify that A^ and Xb are marginal variables when 
OiA + OLB — 1. 

The Monte Carlo results are shown in figure 9 for d ~ 2 and r — 1/2. They 
confirm this behaviour, although the convergence to the asymptotic slope is rather 
slow for the majority species. 
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5. Conclusion 

We have shown, through scaUng considerations and Monte Carlo simulations, that 
when the reaction rate of the single-species diffusion-annihilation process decays as 
t~^, there is a critical value lUc = 1 — (i/2 for d < dc separating reaction- limited 
behaviour for u) > coc from diffusion- limited behaviour for u) < cUc- In the reaction- 
limited regime, mean-field theory is always valid and leads to an w-dependent decay 
of the particle density whereas in the diffusion-limited regime, due to concentration 
fluctuations, the particle density decays as t"''/^ when d < dc = 2 as for a constant 
reaction rate. 

In the case of coupled reactions where the annihilation of one type of particles 
is catalysed by the other, one obtains a time-dependent effective reaction rate which 
is easy to identify when the process is reaction limited, i.e., when mean- field theory 
applies. 

The main effect of decreasing time-dependent reaction rates is to extend the 
validity of mean-field theory below the usual upper critical dimension and to lead to 
unusual mean-field behaviour. 
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